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The eigenvalues of Wilson loop matrices in SU(N) gauge theories in di- 
mensions 2,3,4 at infinite N are supported on a small arc on the unit circle 
centered at z = 1 for small loops, but expand to the entire unit circle for 
+^ large loops. These two regimes are separated by a large N phase transition 

O | whose universal properties are the same in d — 2, 3 and 4. Hopefully, this 

D large ./V universality could be exploited to bridge traditional perturbation 

theory calculations, valid for small loops, with effective string calculations 
for large loops. A concrete case of such a calculation would obtain analyti- 
cally an estimate of the large N string tension in terms of the perturbative 
> scale A SU{N} . 

On 
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I. Introduction 

These lecture notes provide an elementary introduction to the topic 
described in the abstract. Most analytic results are presented as exercises 
- the solutions of most of these exercises can be found in the papers in the 
J> reference list. These derivations are not reproduced in the notes. Rather 

than describing in detail the numerical results, the focus is on the general 
logic of the numerical tests. Again, details are in the references and actual 
numbers and graphs are not reproduced in the notes. These notes reflect 
my personal viewpoint. 

The entire research topic rests on a paper by Durhuus and Olesen, |T], 
from 1981, who found a large N phase transition in two dimensional YM. 
Much of the work on the crossover was influenced by Blaizot and Nowak 
who pointed out almost two years ago that one might view the onset of 
confinement as an analogue of the onset of Burgers' turbulence [2]. 
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2. Abelian Wilson loop operators 

We work in the framework of Euclidean Field Theory, in R d , where 
d=2,3,4. The focus is on pure gauge theory, that is, there is no matter. 

2.1. Free abelian gauge theory 

We start with U{\) gauge theory, a free, Gaussian, theory. The quantity 
below, Z[J^\, contains all of the information about this theory: 

Z[J,] = J [dA,] e -^f dd * F ^f dd ^ A ^ J [dA ^ e -S[A, ]+ if d « x j»A, 

(2.1) 

All vector indices [i{v) are pairwise contracted. The field strength F^ v is 
defined by 

F lw = d ll A v -d v A ll and F% = F^F^ (2.2) 

2.2. Carrying out the integral over 
2.2.1. Decoupling the modes 

The integral is Gaussian and there is translational invariance, so we can 
decouple the modes by Fourier transforming the fields. 

k 

(2.4) 

We now decompose the vectors A^k) 

A,{k) = ^a L {k) + e^{k)a\{k) (2.5) 

with 

6^ = 0, ^(k)ej i (k) = S i t, i,j = l,..,d-l (2.6) 

The quadratic term in the exponent of the path integral does not depend on 
ai{k). Hence, if the current obeys k^J^k) = or, equivalently d^J^(x) = 0, 
one can introduce another weight factor for ai(k) to make the integral per 
fc-mode finite, without Z[J^\ depending on the details of that weight factor. 
One can even drop the aiik) variables altogether. 
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2.2.2. Current conservation 

If one introduces a J M (x) which is not conserved, that is d^J^(x) / 0, 
the integrand does not have a local extremum around which to expand. The 
first variation of the exponent gives 

= g 2 J v (2.7) 

If the above equation for AAx) had a solution, this would imply that 
d^J^x) = by virtue of F^ u = -F Vfl . 

One could interpret the a l integration as forcing the longitudinal compo- 
nent of to vanish. The theory is simply restricted to probing by conserved 
external currents only. 

More generally, only gauge invariant observables are meaningfully de- 
fined by the theory. The gauge transformation is 

A^k) - A^k) - ik^xik) (2.8) 

It only affects ol; gauge invariant observables do not depend on at,. 
The most localized conserved current one could imagine is 

pi A 7 

J li (x) = J dr5 d {x-z{r))^ (2.9) 

where z^(t) is a closed contour described in R d by r varying from r = to 
t = I. The current is defined by the curve itself, not by its parameterization. 
We fix the latter by 



dV 2 



dr 

This makes / the perimeter of the curve, since 



(2.10) 



/' 

Jo 



The curve had to be closed to ensure current conservation: 

d, J M (x) = jj dr^^ix - z{t)) = - jj dr£5 d (x - z{r)) = (2.12) 

5 d {x - z(0)) - 5 d (x - z(l)) = (2.13) 

In Fourier space 

J (k) = ( d d xe ikx JJx) = f dr^e ikz{T) (2.14) 
J Jo dr 

Current conservation now is seen from 

kMk) = -i [ l dr^e lkz ^ = (2.15) 

JO UT 
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2.2.3. Wilson loop operator 

The Wilson loop operator is W[A] = e i f ddxJ n A n = j § dz^A^z) ^ ma ki ng 
the independence on curve parameterization explicit. Its average is 

(W) = (2.16) 

2.2.4. Overcoming basic problems 

We want to calculate (W). We face some problems: 

• There is no weight for A^(0). Luckily, J^(0) = so J M is not coupled 
to this degree of freedom and we can forget about it, taking the inte- 
gration over it to cancel between Z[J^\ and Z[0]. Wilson loops appear 
to be "infrared safe" . 

• The integral is also over ai(k) and is unbounded in that direction; 
as mentioned, we can fix this by putting in a weight; this is done by 
multiplying the integrand by 

e -^/^^) 2 = -^iS 0tk*a L (k)a L (-k) (217) 

The ao dependent terms cancel between Z[J^\ and Z[0] leaving (W) 
ao independent. 

• There is an infinite number of k- values, with ~ amplification as 
k 2 becomes large. The factorized path integral consists of an infinite 
number of finite factors, and the product over all of them does not 
converge. This is solved by introducing an ultraviolet cutoff A and 
restricting the integration variables by 

k 2 < A 2 (2.18) 

Consequentially, the averaging process over the A^kys is now only 
sensitive to J^(k) with k 2 < A 2 . The path integral does not affect the 
dependence of W on J^(k) with k 2 > A 2 . 

• Ju(x) is a distribution: it is zero for x ^ z{t) Vt and infinite if for 
some r we have x = z(t). A Gaussian integral is bound to give (taking 
into account translational invariance) a Gaussian answer: 

e -| / d d xd d yj \ 1 (x)J v (y)G y. v (x-y) ^ 19) 

Unless G^ u (x — y) makes this irrelevant, we have to confront products 
J fM (x)J u (y) at coinciding points x and y. There is no general way to 
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give meaning to the product of distributions at the same point. This 
problem is now naturally solved by simply changing our definition of 
the current, first in Fourier space, by setting 

Jn(k) = for k 2 > A 2 (2.20) 

The new current is still conserved: 

4(x) = [ f\,e- ikx ( dr^e 1 ^ (2.21) 
* v ' J k 2 <A 2 (27r) d Jo dr y ' 

J A (x) is no longer sharply localized at the curve; rather it is just 
peaked at it and spread out over a distance 5 ~ ^ in its vicinity. 
Actually, because of my simple choice of a sharp momentum cutoff 
this localization is not good enough for all cases, as the decay away 
from the curve is slow and there is an oscillatory behavior of high 
frequency. 

In order to distinguish the curve's shape from that of an amorphous 
blob, we need Al S> 1. In principle, we want to rid ourselves of the 
dependence on the large number Al by taking it to infinity. Whatever 
we can make sense of in that limit is a universal feature of the theory, 
since many different schemes of overcoming the problems listed would 
produce identical results when the cutoff is removed. We now turn to 
see what we can do; this depends on the dimension d. 



2.3. Circular loop 

Exercises: 

Consider a loop given by a circle of radius R = 1/(2tt) in the 1-2 plane. 
Take d = 4. 

1. Calculate (W) for J A . Answer has the form (W) = e~ g2l l 2 with 

/ ~ -(AR) 2 (\\og^J 2 {RA^l) (2.22) 
Jo 

where J n are Bessel functions of integral order n. 

2. What is the leading behavior of / as AR tends to infinity ? Answer: 

/ = c (AR) + lower orders (2.23) 

For obvious reasons the first term is known as the "perimeter diver- 
gence". From the calculation it becomes obvious that for a curve of 
different shape, but same perimeter, this divergent term would be the 
same as for the circle. 
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3. Consider now a change in the definition of the current J^; instead 



j A. 

of (|2.20b we use a j£' with A' < A. We still have MR > 1, say 
MR = 100, but A > A' and we now take AR to infinity at MR fixed 
at some number much larger than one. What happens now ? Answer: 
The perimeter divergence disappears, and gets replaced by a finite 
perimeter term given by MR times some number as the leading term 
in an expansion in ^t^- 

Take now d = 3 and consider the same questions. The main difference is 
that now the perimeter divergence is just logarithmic. Finally, take d = 2. 
Now there are no divergences. 

The final conclusion is that the best place to start learning about Wilson 
loops is in two dimensions. There one can consider true, infinitely thin, 
curves. It may be possible to study loops also in higher dimensions, but 
this may require to "fatten" the curves by a finite amount, of relative order 
Ml. 



2.4- Geometric meaning ofW 

2.4.1. Matter 

In full QED, photons interact with Dirac fermions of charge g. The path 
integration is extended to include the Grassmann variables ip a {x) 1 ip a {x) and 
the action iSf-AJ is changed by an additive term 

S+ = J d d x [tP(x)^D^{x) + m${x)1){x)] ; D„ = d*- iA^x) (2.24) 

[The Dirac indices a, (5 are silent and summed over. They will not appear 
again.] This coupling preserves gauge invariance with being replaced by 

vvv = 4 ( 2 - 25 ) 

Now, d^J^ is not zero for an arbitrary ip,ip, but, it is zero when the fields 
ip,ip both satisfy the extremum condition along with An, and an expansion 
of the entire path integral in g can be defined. If one instead integrated 
out ip,ip first, the result sf^lA^] is a gauge invariant addition to SLA^]. If 
one instead integrated out A^ first, the component of that integration 
effectively puts a constraint on the remaining ip, ip that enforces = 0. If 
one has added a weight factor for a^, for any positive oo no constraint on ip, ip 
gets generated and only a subset among all possible A^, tp, ip observables are 
deemed physical. That is the subset of gauge invariant observables; they 
would be independent of the parameter ao- 
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2.4.2. Holonomy 

For g = 0, = and we know that generates translations: 

e a ^ip(x) = %j){x + a) (2.26) 
In particular, for a closed curve we have 

e§ ig ^{x) = ij){x) (2.27) 
Replacing by in the above equation gives 

J dz^-iA^z)]^ = w *^ x ) (2.28) 

So, when ip(x) is "transported" round a closed curve it accumulates an 
additional phase factor (holonomy), given by the Wilson loop operator as- 
sociated with the curve. 

2.4.3. Infinitesimal loops 

For an infinitesimal loop, bounding a flat surface element 5a one has 

W « 1 + iSa^F^ (2.29) 

The role of S[Afj] is seen to provide a bound on the fluctuation in W when 
the loop is very small. This works well in d = 2, but, as we have seen, has 
problems in d > 3. 

2.4.4. Some dimensional analysis 

For any d, W is dimensionless. This implies that has dimension 
1 (this is mass dimension, the inverse of length) and therefore F^ u has 
dimension 2. S'f^l^] is in the exponent, so also needs to be dimensionless. 
Therefore, the dimension of the coupling g 2 is 4 — d. 

After averaging, one gets for a small loop, schematically, 

1 - (W) ~ (5aF) 2 (2.30) 

Here we ignored the perimeter term, although it is divergent for d > 2. 
Let 51 be the linear scale of the loop. The gauge action gives 

^{5l) d F 2 ~ 1 (2.31) 

Hence, with 5a ~ 5l 2 , 



1 - (W) ~ (<5/ 2 F) 2 ~ ^((M) 4 ^ 



(2.32) 
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For a small loop this would go to zero as the loop shrinks, except at d = 4, 
where the theory is scale invariant in the present approximation. Once the 
approximation is improved upon, scale invariance is lost even in d = 4 and 
g 2 starts depending on 51, logarithmically for small 51. If the gauge theory 
is nonabelian, one has g 2 (5l) ~ — 1 / log(Ag U ^51) for A SU ( N ^5l — > and 
one recovers (W) — > for 5/^0 even at d = 4. 

However, for d > 2, the above effect is masked by the perimeter diver- 
gence. 

2. 5. Summary of section 

• Wilson loops associate a phase factor with parallel transport round a 
closed curve 

• The gauge action is defined so as to suppress fluctuations of phase 
factors associated with infinitesimal loops away from identity 

• For small loops the phase factor is close to identity, but this is masked 
by a perimeter divergence for d > 2 

• For d = 4 to behave similarly to d = 2,3 we better focus on non- 
abelian gauge theories, whose effective coupling tends to zero at short 
distances 

3. The nonabelian holonomy 

3.1. Definition 

The nonabelian gauge group has a Lie algebra Q of real dimension n 
and generators labeled by i = l,..,n. We restrict our attention to Q = 
su{N) with n = N 2 — 1. There are now n vector fields A^(x). su(N) is 
compact and has a discrete infinity of irreducible finite dimensional, unitary 
representations, in which the generators are represented by linear unitary 
operators acting on finite dimensional Hilbert spaces. Each representation 
is labeled by R, the dimension of the associated Hilbert space by dR and the 
operators representing the generators by Choosing a basis, makes the 

generators traceless hermitian matrices with entries denoted by , a,b = 
1, .., dR. We adopt the convention 

TrT (R)i T (R)j ^ 5 ij ( 3 .!) 
so that the structure constants Cijk are totally antisymmetric: 



[J>(.R)^ rp(R)j^ _ i( jijkrp(R)k 



(3.2) 
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For any R, we define a covariant derivative acting on a matter field 
ipl R \x) both on x and on b: 

[D^] a (x) = [d*5 ab - iT^Afa^Mx) (3-3) 

Parallel transport round a closed curve has a holonomy given by the 
action of a unitary matrix, representing a group element W in the repre- 
sentation R: 

W R {x) = Pe i $ ds M*)* It >' (3.4) 

3.2. Gauge invariant content of holonomy 

The symbol P indicates that the product is ordered round the closed 
path, starting at x and ending at the same point x. Ordering is necessary 
as the matrices in the exponent, at different points on the path, do not 
commute with each other. A finite gauge transformation acts by an x- 
dependent group element g(x) 

W R (x) - g^(x)W R (x)g^(x) (3.5) 

where g( R \x) is the unitary matrix representing g(x) in the irreducible 
representation R. Hence, the gauge invariant content of the holonomy is 
contained in the infinite collection of numbers 

X r(W) = Tr W R (x) (3.6) 

The notation indicates that the dependence on x disappears, and these 
are indeed numbers associated with the closed curve, without selecting any 
particular point on it. W is the abstract group element associated with 
parallel transport round the curve from x to x. 

Viewed as functions on group space, the Xr(W) are a basis of the space 
of all class functions on the group. Any class function is a smooth function 
of the eigenvalues of 

W f (x) = Pe i $ dz » A ^ T(f)j (3.7) 

where / denotes the fundamental representation; for su(N), it is by hermi- 
tian, traceless N x N matrices. The eigenvalues are TV" points on the unit 
circle e* 7a , 7 a real, constrained by Il^o 1 e * 7a = 1- Wf is an TV x N unitary 
matrix of unit determinant. 
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3.3. A probabilistic view 

As is well known, one can write a gauge invariant Lagrangian for the 
fields A l ^(x), but it is nonlinear. For a small loop one has 

W R ~ 1 + ida^F^TW" (3.8) 

where 

= d,A k u - d u A* + C^A^Ai (3.9) 

The action S[A l ] still consists just of a suppression factor in the traces of 
holonomies round small loops: 

S W = 4^2 / ^« ( 3 - 10 ) 

The integrand is a class function, so the suppression is on the eigenvalues 
of the small loop Wf. The action tries to make Wf close to the N x N 
unit matrix for small loops. In other words, the action tries to make the 
holonomy of tiny loops close to the identity in the group. 

One could imagine integrating out in the path integral all A 1 (x) with 
the constraint that, for a chosen fixed loop, its holonomy W is kept fixed. 
This would produce a positive function (we are setting to zero a potential 
second term in the action, that is possible at d = 4, and which is imaginary) 
since the integrand is positive. Normalizing, we would get a probability 
distribution on the group, P(W). Gauge invariance tells us that P(W) 
must be a class function and that the measure of integration on W must be 
the Haar measure. P(W) contains all the information needed to compute 
all averages of the form 

(X R1 (W) X R 2 (W)....) (3.11) 

Since P(W) can be linearly expanded in xr(W) the coefficients in that 
expansion hold all the information determining the above moments. 

One can think about S^-AJJ as defining a highly peaked, "bare" distribu- 
tion Pq(W) for every tiny loop; the peaking is controlled by g 2 . The weight 
generated by the Pq(W) for all tiny loops W collectively produces the above 
P(W) for a macroscopic loop. 

This picture has ignored the problems we have encountered in the case 
where the group is U(l). These problems have not gone away, since for 
g 2 — > we have just iV 2 — 1 non-interacting photons. However, there are no 
problems at d = 2. 
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34- S[Al} from P (W) 

Let us start with Pq(W), from which we wish to derive S[A l ] and P(W). 
To make an action out of Pq(W) one needs to make sure that the collec- 
tion of infinitesimal loops that is selected includes every A 1 (x) (all x,/j,,i). 
There is much freedom in the choice of Pq{W) itself, as only the behavior 
in the vicinity of W = 1 matters. A very natural choice is provided by the 
heat-kernel function: One writes down a diffusion equation for W which is 
consistent with the homogeneity of the group manifold: 

§^Po(W;t) oc V 2 w P (W;t) (3.12) 

t > is dimensionless and is the Laplacian on the group manifold: It 
is defined by being the ordinary Laplacian in the tangent space at W = 1 
extended by conjugation to the entire group. maps class functions into 
class functions. The diffusion constant, which could be absorbed in t, is 
a number chosen by some convention, to make the formula for Po(W,t), 
below, correct. The initial condition for the diffusion equation is taken as 

P (W;0)=5 U ^(W,1) (3.13) 

Thus, Pq(W; t) will be a class function for all t > 0. It is explicitly given by 

P (W;t) = Y,dRXR(W)z-^ C2{R) (3-14) 
R 

Here the sum over R is over all irreducible representations of SU(N) and the 
number C2(R) is the value of the quadratic Casimir operator on the repre- 
sentation R in a specific convention for the normalization of the generators. 
For all t > and W we can view Po(W;t) as a probability density since 
Pq(W; t) > and the diffusion equation preserves the initial normalization. 

The group elements W ~ 1 are identified by Wf = e tH ~ 1 + iH, where 
H is iV x N, hermitian and traceless. For t< 1, the matrix norm of H is 
small with high probability. We now see that we can identify 

H ~ Sa^F^T^ (3.15) 

where the small loop in question is the boundary of the little surface element 
5<Jfiv We ensure that all Aj^(x) participate by including all rotated and 
translated copies of da^u, which amounts to multiplying ^(W 7 ,*) factors 
over all \i > v and x. 

Taking t « 1 we can use the tangent space approximation to the 
diffusion equation and we get the standard YM action. 
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3. 5. Summary of section 

The holonomy is a geometric construct and is the central object of non- 
abelian gauge theory. Both the action and all physical observables are given 
by holonomies. 

One can view the basic problem of Euclidean non-abelian gauge theory in 
R d as the calculation of the probability distributions of holonomies for arbi- 
trary loops given the probability distribution of a complete set of holonomies 
of infinitesimal loops. By a complete set, one means a set of loops whose 
holonomies determine the connection A 1 {x), up to gauge transformations. 
In turn, the A^{x) determine all holonomies and are the variables of inte- 
gration in the path integral. 



4. Two dimensions 

In two dimensions the calculation of P(W) for a macroscopic, non- 
selfintersecting loop, given Pq(W), is straightforward. One needs to tile 
the area inside the loop by an exactly fitting cover made out of elementary 
microscopic loops. So longs as such a cover exists, the shape of the loop 
does not matter; just the total enclosed area affects P(W). The result is 

P(W) = P (W;r) (4.1) 

where r is the area A enclosed by the loop in units of A = g 2 N: 

r = XA (l + (4.2) 

Therefore, as we scale the loop up, its holonomy W spreads from the vicinity 
of unity for a small loop, to visiting the entire group evenly at r = oo. 



4-1. The Durhuus-Olesen JI]/ non-analyticity at N = oo 

I shall refer to this non-analyticity as a "phase transition" . This phase 
transition occurs at ./V — > oo at fixed A.4 when r increases through the value 

T C = 4. 

There are many ways to see the phase transition. I choose one of the 
simplest, employing the average characteristic polynomial of Wf. 



4.1.1. Average characteristic polynomial [3 



Exercise: Equation (4.1) gives for a non-selfintersecting loop enclosing 



an area equal to r in dimensionless units: 

(XR(W(r))) = due-™*™ 



(4.3) 
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Exercise: Define 

(z,t) = (det(z - W/(t))) (4.4) 

Show that the above polynomial in z generates all the {xr(W(t))) with 
totally antisymmetric R (single column Young pattern). 
Exercise: Define 

6W(z r)= 1 -- 1 Z u 5) 

Prove that 

^ N \y,r) = ^ N \-ey,r) (4.6) 

with real y, satisfies: 
with initial condition 

^ N \y^) = -\t^ y - (4.8) 



Exercise: Show that the Burgers' equation in (4.7), with the above initial 
condition, produces, in the limit N = oo, a shock at y = when r reaches 
the value r = 4. For r < 4 the solution is smooth in y. 

We see that ^ plays the role of viscosity in the Burgers equation. The 
zero viscosity limit is singular, producing a "breaking wave" at r = 4. The 
singularity is absent at any finite N; hence the term "infinite N phase 
transition" . 

Note that the critical value r = 4 is determined by the initial condition 
and not by the equation. Also note that deforming the initial condition 
would not eliminate the appearance of the shock, except when the initial 
condition is changed in a drastic manner. So, while the appearance of a 
shock is a likely event, the exact value of r where it happens varies from 
case to case. 

Exercise: Explain Figure [TJ specifically, what the arrows mean and why 
the caption is justified. 

4.1.2. Motion of zeros |4j 

As a monic polynomial, (—l) N ip( N \z,T) can be parametrized by its 
zeros. 

Exercise: Prove that all the zeros z a (r), a = 0, .., N — 1 are on the unit 
circle. 
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Shock Formation 





Fig. 1. The making of a shock. 



Exercise: Write z a {r) = e lda ^ and work out the equations of motion 
for the angles 9 a (r): 



d9 a 1 v-^ ,0 a — Ob 



dr 2N 



b, bj^a 



The initial condition is at a singular point: 

0a(O) = 



(4.9) 



(4.10) 



Thus, the partial differential equation (4.7) with our specific initial condition 



can be reduced to a system of N first order differential equations. 

Using these equations one can evaluate the motion of the zeros for rCl 
and for r»l. 

Exercise: Make the ansatz 



(4.11) 



9 a (r) = 2i M /- 



and solve the equations of motion for the 6 a {r) to leading order in r <C 1. 
The result is that the rj a are the zeros of the Hermite polynomial of order 
N: 

H N ( Va ) = 0, a = 0,l,...,N -1 (4.12) 



Exercise: Prove that for all r > 

N-i 



(4.13) 



a=0 
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and that y pairs of angles have the same absolute values and opposite 
signs. Also, show that if N is odd there is one angle that stays at for all 

T. 

Exercise: Find the asymptotic behavior of the zeros as r — > oo. The 
result is 

a {r = oo) = ^ (a - ^) e a (4.14) 

At infinite r the zeros are uniformly spread round the unit circle. 

Exercise: How do the zeros approach their r = oo values ? Set # a (7~) = 
@a + 50 o (t). The result, to leading order, is 



S0 a ( T ) ~ -2e-^ (A '- 1) sine a 



(4.15) 



AT 3iV 
4 ' 4 



- — . — move the fastest at large r. 



Note that the zeros with a 

Exercise: Let N S> 1 and fixed. Calculate the asymptotic behavior of 
the pair of zeros closest to —1 on the unit circle at r = 4. The answer is 



ZM 



exp 



± 



3.7*' 



(4.16) 



The number 3.7 is an approximation. 

Exercise: Let N ^> 1 and fixed. Let | = 1 + -^7. Find f such that 
zm(t) has a finite nontrivial dependence on a of the form 



zm ~ - exp 



/(a)i 



3 



(4.17) 



as — > 00, with /(0) ~ 3.7. The answer is = 1/2. 

We see that the motion of the extremal zeros becomes nontrivial in the 
limit N — > 00 if we express them in terms of specific scaling variables. The 
associated exponents of iV are 1/2 and 3/4. 



4-2. Zeros and eigenvalues 

At large enough N, the set of zeros z a {r) is a good approximation to the 
set of most likely values of the eigenvalues of the fluctuating matrix Wf. 
Too see this we need to compute the single eigenvalue density pn(6', W). 



Pn(9;W) 



1 

N 



(4.18) 



where the the eigenvalues of an instant of Wf are e* 7a , a = 0, .., N — 1. 

Exercise: Compute pn(9,t). Start by expanding det(l + uWf) / det(l — 
vWf) in characters, then take the average and next study the limit u — > — u. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Fig. 2. Plots of the density pn(9) (oscillatory curve) together with the positions of 
the angles of the zeros 8 a (vertical lines) for r < 4 (left) and r > 4 (right), N = 10 
(top), and N = 50 (bottom). 

At N — oo, as first shown by Durhuus and. Olcscn pQ, Poo 

(9, t) has a gap 

centered at 9 = ±7r for r < 4, which closes for r > 4. The exact formula 
for finite N shows that there are N oscillations modulating the N = oo 
form and, when the N = oo case has a gap, the finite N density has a 
non-vanishing tail, exponentially suppressed with N in the gap region. 

The peaks of the oscillations can be interpreted as the most likely lo- 
cations for the fluctuating eigenvalues of Mt, with phases j a . The plots in 
Figure [2] show how they compare to the locations of the zeros z a . We may 
loosely refer to the zeros z a 's as the eigenvalues of Wf or even W. 

The behavior of pn(9,t) close to critical r and for 9 close to ±7r is 
universal at large N; this statement can be made precise by defining "close" 
for the two parameters 9 and r with the help of two scaling variables £, a and 
the associated critical exponents of N, 3/4, 1/2 respectively. The extremal 
zeros zm reside in the critical regime of 9 when r is also in the critical 
regime. 

4-3. Summary of section 

For a small loop, the eigenvalues of Wf are all concentrated around 
unity. As the loop expands the eigenvalues expand round the unit circle 
until they cover it completely. Very large Wilson loops have eigenvalues 
randomly distribute round the unit circle. There is a crossover between the 
two regimes. This crossover can be seen in various ways: The evolution of 
the average characteristic polynomial as a function of y, where z = —e y with 
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y-real, shows a sharp behavior at y = when the size of the loop is in the 
crossover regime. The extremal eigenvalues close the gap in the crossover 
regime. When taking N — ► oo the crossover turns into a phase transition, 
and for large enough N there is a regime of universal behavior. 

The main qualitative change between the perturbative regime for small 
loops and the non-perturbative regime for large loops is that the eigenvalues 
cover just a small arc round z = 1 in the former and the entire unit circle 
in the latter. The crossover is when the entire unit circle first gets covered. 
The zeros of the average characteristic polynomial go as e t ^ rria for small r 
and like exp {i [@ a — 2e~ aT sin © a ]} for large r, where a is the fundamental 
string tension and the G a cover the circle uniformly. 

Our main hypothesis is that this behavior also holds in d = 3,4 once 
holonomies are properly defined, eliminating the perimeter divergences. If 
this is true, one can hope that in d = 4 the regimes of confinement and 
perturbation theory are separated at large TV by a narrow regime in which 
one has universal behavior. Further, I hope that this will allow to actually 
match perturbative to non-perturbative behavior in d = 4 leading to a way 
to estimate the string tension in four dimensions in terms of a perturbatively 
defined scale, h-su(N)- 

In short, the holonomy evolves with scale from exploring an arc round 
z = 1 to filling the entire unit circle and the crossover between the two 
realms narrows as N — > oo. For iV S> 1 the crossover can be described 
universality, using a simple random matrix theory effective description of 
W f . 

5. Higher dimensions 

In d > 2 our hypothesis needs to be tested numerically. I shall set up 
first the general (that is, for d = 3 and d = 4) strategy for carrying out a 
test. 

This strategy has been implemented fully in d = 3 and partially in 
d = 4. d = 3 is cheaper in computer time than d = 4 and there is one 
additional simplification: there are no "corner divergences" . We have seen 
that at d > 2 there is a perimeter divergence. Only the existence of a 
piecewise continuous tangent to the curve enters in determining the term 
of highest degree of divergence. If the curve has kinks, that is, it's tangent 
has isolated discontinuities, sub-leading logarithmic divergences localized 
at the kinks develop. At d = 3 the perimeter divergence is logarithmic, so 
kinks produce no sub-leading divergences. For numerical tests one employs 
typically hyper-cubic lattices and closed contours on the lattice must have 
kinks. In d = 2 the curve only needs to be continuous, even the existence 
of a tangent is not required. 
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5.1. Smearing fjjj/ 

We need to eliminate the perimeter divergence. As anticipated, we do 
this by "fattening" the loop; this has a major effect on small loops and a 
limited effect on large ones. We need to extend the ideas from the abelian 
case to the non-abelian one. We introduce an extra parameter p > [do 
not confuse the parameter p with the eigenvalue density pn(0;W)\, and 
make the A 3 (x) variables, that enter the definition of the Wilson loop, p 
dependent. In this sense, one may think about p as adding a fifth dimension, 
taking R d to R d x R + . But, it is only the Wilson loop that feels p. The 
pure Yang- Mills theory being probed lives still at p = 0. We define new 
variables A 3 {x,p) by a diffusion-like equation in p: 

d ^j^ = [D fl F, u (x,p)Y (5.1) 

where is the covariant derivative in the adjoint representation with re- 
spect to A J Jx, p) and 

F^(x, p) = d^Al(x, p) - dvAfa, p) + C^A%x, p)A k u (x, p) (5.2) 
The equation is supplemented by the initial condition 

^(x,0) = 4(x) (5.3) 

where the A^x) are the fields in Sf-A^J and are the integration variables of 
the path integral. 

The smeared Wilson loop is defined by 

W R = PexpUdZpAfatptTW'] (5.4) 

Although we do not do this here, it is appealing to think of replacing the 
constant p above by a function p(r), where the loop is described by z^{t). 
It is important to realize that the behavior under gauge transformations in 
R d is preserved. 

How does this fatten the loop ? At leading order in perturbation theory 
we just have iV 2 — 1 photons and the equations linearize. In Fourier space 
we have 



dp ~ l 



(5.5) 

up 

with 



AUk,0) = AUk) (5.6) 
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The longitudinal component of A^k, p) does not vary with p but it does no 
appear in the action and does not couple to the Wilson loop either. Each 
of the 3 transverse components, for each j = 1, N 2 — 1, goes as 

a{(k,p) = e- pk2 a{(k) (5.7) 

Hence, the coupling of the momentum modes of A^{x) with k 2 p S> 1 to the 
closed curve is exponentially weakened. We have not literally fattened the 
curve - it stays a mathematical one dimensional curve, but the path integral 
integration variables do not resolve the curve well over distance smaller than 

To avoid losing the dependence on shape for small curves altogether, we 
can make p depend on the scale of the loop, perhaps most fittingly on its 
perimeter (assuming, say, that the loops is an arbitrary, but smooth and 
small, deformation of a circle). For large curves there is no danger of losing 
sight of the loop shape so we imagine choosing 

72 

(5.8) 



r [lA'] 2 + c 

where c is some large number, say c = 20. 

The above construction is quite ad-hoc, but we hope that it will allow us 
to see the universal large N crossover we expect in analogy with the d = 2 
case and that the latter will not depend on the details of our "fattening" 
procedure. 

The main advantage of the "fattening" procedure we chose is that it can 
be easily extended to lattice gauge theory, so we can test our hypothesis 



outside perturbation theory. The point is that the right hand side of (5.5) 
is the variation of the classical action, and therefore has a natural lattice 
counterpart, where the continuum integration variables A^ are replaced by 
the lattice integration variables U^{x) [a SU{N) N x N matrix attached 
to the link connecting the grid vertex at x to its nearest neighbor in the p 
direction] and the continuum action by the lattice action. The lattice action 
itself is literally constructed as a product over all identical, sharply peaked, 
probabilities for the holonomies going round all the most elementary loops 
on the hypercubic lattice [going round small two dimensional flat squares 
called plaquettes], so the basic concept on the lattice is the same as in the 
continuum. In practice, the differential equation in p is discretized in the 
lattice context because it has to be solved numerically, being nonlinear even 
in the continuum; there are various variants of doing this and the precise 
procedures are known in Lattice Gauge Theory as "smearing" methods. 

Note that the smearing method will also eliminate "corner divergences" 
so it is adequate also for d = 4 dimensions, not just d = 3. 
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5.2. The lattice test 

The lattice is taken as a <i-dimensional hypercubic grid of finite extend 
in each direction and connected as a torus in each direction. On each lattice 
link one has the link variable U^(x) described above. The links are all given 
an orientation and parallel transport is now taking a matter field from one 
site to the neighboring one by multiplication by Unix) if the link is traversed 
in the positive direction, or U^x) if the link is traversed in the negative 
direction. The collection of all U^x), Vx, fx can be thought of as defining the 
d operators e Dfl . Transporting from site to site one can parallel transport 
along any path made out of consecutive links. The parallel transporters 
round the elementary plaquettes, identified by a site x and directions fj,, u 
are denoted U p where p identifies the plaquette. The lattice action is 

S[[/ M (z)]oc^Tr[l-^(x)] (5.9) 

v 

One uses the notation (3 for 4^ and b for 

Given a set {U^x)}, one produces a set of smeared SU(N) matrices 
{U^(x,p)} from which one constructs smeared Wilson loop matrices Wf 
using the smeared link matrices as elementary parallel transporters. For 
definiteness, we shall only consider square loops, in a lattice plane. These 
loops have a side consisting of L links, where L is an integer. 

There are no dimensional parameters on the lattice; dimensional analy- 
sis in the continuum is just an aid in carrying out various scaling transfor- 
mations. There is no intrinsic meaning to dimensional quantities in Field 
Theory. 

The total number of sites of the lattice is denoted by the integer V and 
is referred to as the "volume". One needs to make V so large that the 
results obtained can be extrapolated to the limit V — > oo. This is referred 
to as the thermodynamic limit and, excepting computational cost, causes 
few problems of principle in pure non-abelian gauge theories because they 
have a mass gap. 

On the lattice, everything that was said before, which did not quite 
hold because of various ultraviolet problems, now holds in a precise sense. 
We do have a well defined problem in classical statistical mechanics. The 
continuum results are obtained by taking the continuum limit. One must 
make sure that one only poses questions that have finite non-trivial answers 
in the continuum. The ideology behind this is the concept of Field Theory 
universality, which says that physical continuum answers will not depend on 
any details of the lattice definitions and will admit an asymptotic expansion 
that coincides with the asymptotic expansion in the continuum, defined by 
renormalized perturbation theory via Feynman diagrams. 

Exercise: Show that the above works in d = 2. Smearing is not needed. 
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5.2.1. The ingredients of a lattice simulation 

The heart of the process is a code, a computer program which imple- 
ments a Monte Carlo process that generates sets {U^x)}, "configurations", 
with probability 

e -s[u,( x)] (510) 

The basic idea behind these codes is that they implement a Markov chain 
made out of simple steps to which one inputs a stream of pseudo-random 
numbers generated in one of the better ways known. It is assumed that 
the numerical errors induced by the stream being pseudo-random rather 
than truly random are far below the statistical errors induced by sampling. 
Sampling is used to estimate averages of the interesting observables, Wilson 
loops in our case. One is restricted to finite samples and one wants to get 
good estimates at minimal cost of computation time. 

Thus, one has a computational method of finding out the numerical 
value of the average (xR^Wf)) f° r an Y lattice loop [in our case, square loops 
of side L] to some accuracy. 

There are several parameters one can set: N, the group size; R, the 
representation; b, the lattice coupling; L, the loop size and p, the smearing 
parameter. There are more parameters one is free to set in the code, that 
impact efficiency, controlling several aspects of the algorithms and the pre- 
cise form of the Markov chain. We shall assume they have been fixed and 
are not changed during the simulation. 

The result of the simulation would then be a table of numbers for 
(xr(W)) one entry for each set (N, R,b,V, L, p). If one intends to also 
test for Field Theoretical universality, one can vary over different P$(W) 
functions, changing the action S\U^{x)\. 

5.2.2. Analysis 

Once the tables of numbers have been produced (the "data" ) one needs 
to analyze this data. One part of the analysis is purely statistical and 
determines the reliability of the data, expressed in one of the standard 
ways one quantifies statistical uncertainty. Loosely speaking, this stage 
determines the "errors" . 

Here I shall only focus on the second stage of the analysis, the extraction 
of estimates in the continuum limit. Since this stage may involve the data 
in very nonlinear ways, another step of statistical analysis is needed, to 
determine how the errors propagate to the final answers. 

I shall not describe any of the error analysis, nor the algorithm used to 
generate the data. 
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5.2.3. How our hypothesis is tested |6J 

To get to the continuum limit one needs to take b — > oo and V, L — ► oo 
at fixed b in a correlated manner. One tries to make V large enough from 
the start, so that one can forget about it. This is easier as ./V gets larger 
because of a phenomenon known as "large N reduction" . Large N reduction 
implies that increasing N reduces the size of finite volume corrections. 

To determine how the limits L, b — * oo have to be correlated one keeps 
a selected physical quantity fixed. For example, one could calculate the 
lattice string tension £(&, N) 

E(6,JV) = -lim bg( ^ (L)) (5.11) 

One then defines a scale, a dimensionless number (dropping the explicit 
mention of the dependence on N) A (b), by 



A(6) = y/E(b,N) (5.12) 

It is very common to assign a dimensional number to A (6) called the "lattice 
spacing" a(b) by 

where I just picked some approximate value for the QCD string tension as 
(440 MeV) 2 . For large L 



log{W f (L)) ~ -S(6,A^)L 2 (5.14) 

and we write 

S(6,iV)L 2 = [La(b] 



A(b) 



o(6) 



(5.15) 



Both factors on the right hand side are physical dimensional quantities in 
the continuum limit. 

It is a fact that as b — > oo, N) goes to zero, so a(b) goes to zero 
as b — > oo. This fact is relatively easy to understand. What is less trivial 



to establish is that the left hand side of eq. (5.11) is nonzero. The latter 
fact says that there is confinement on the lattice. To actually establish 
confinement in continuum one needs to also show that a(b) goes to zero in 
the precise way predicted by continuum perturbation theory. 

With a(b) established, the continuum limit is defined by taking L — > oo 
in such a manner that I, defined by 



I = La{b) 



(5.16) 
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is kept finite at a physical distance. That the distance is "physical" is usually 
expressed by converting A(6) to a(b) and attaching units to I. I shall use 
this language from now on, but, remember, the entire calculation has no 
idea what an MeV is. 

For each b we compute, with L = l/a(b), at I held fixed, the eigenvalue 
density 

p N (0;W f (L,p(L(b),b)) (5.17) 

Here 

I 2 L 2 
6) = a 2 (b)[l 2 + ca 2 (b)\ = a 2 {b)[L 2 + c] (5 ' 18) 

where c is a pure number, say 20, that is, significantly larger than 1. 

In principle, we would like now to carry out the following ordered set of 
steps: 

1. Select an I and a N. 

2. Select a sequence of increasing b values, for which V is large enough to 
be assumed infinite, paired with a sequence of L's such that eq. (5.16) 
is obeyed. 

3. For each pair (b,L(b)) in the above sequence calculate the eigenvalue 
density p N (9;W f [L(b), p(L(b),b))) 

4. Extrapolate b — > oo and determine the continuum eigenvalue density 
Pn(0,1). 

5. Repeat the steps above, still keeping I fixed, but increasing N 

6. Obtain, by extrapolation to N = oo, p oc (6,l). 

7. Repeat the above steps, now varying I. The range of I one should use 
should include the transition point l c 

The hypothesis will be supported by the test if we find a finite l c equal 
to c GeV" 1 , where is c is a pure number of order 1, say between 0.1 and 
10. The part of the hypothesis that is being tested is just the existence of 
the transition, defined as the demarcation point separating a Poc(&,l) with 
a gap around 6 = ir for I < l c from a Poo(8,l) with no gap at 9 = tt for 
I > l c . 

Next we need to devise a test of large N universality. This is done as 
follows: 

We start by defining the lattice average characteristic polynomial: 



N {b,L,p(L,b)) = (det (e% + e-%W f (L))) (5.19) 
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We would like to take b — > oo. We already know from the first part of 
the hypothesis test where the phase transition is expected to occur on the 
lattice, and have checked that the critical size of the smeared loop has a 
reasonable continuum limit. 

We are interested in the region y ~ 0, N — ► oo and I ~ l c . We now 
proceed by ordering the ./V — > oo and b — ► oo limits differently This amounts 
to making an extra assumption: We assume that the infinite N limit of the 
continuum limit is the same as the infinite b limit of the infinite TV" limit of 
the lattice theory. So, at fixed b (which fixed L and p), we take N to infinity 
first. The successful result of the first part of the test is the main reason 
why we are willing to add this assumption. The reversal of limits simplifies 
the procedure significantly. 

To test for the universality component of the hypothesis we need to 
identify at fixed L and at infinite N a critical coupling b c (L) where the 
spectrum of Wf(L) just closes its spectral gap at eigenvalues equal to -1. 
a{b) is monotonically decreasing as b increases. So, the physical size of 
the loop I = La{b) shrinks as b increases. Dilating the loop corresponds 
to decreasing b. Varying b at fixed L will take us through the transition. 
Looking at plots of pn{0) for Wf, one can observe this. A more quantitative 
method is as follows: 

Take the numbers Ojv(6, L, p(L, b)) and expand in y around y = 0. I 
suppress the dependence on L, as it is held fixed. 

N (y, b) = C (b, N) + d(b, N)y 2 + C 2 (b, N)y 4 + ... (5.20) 

Define 



n(b,N) = " uyv £;"^ ■ ( 5 - 21 ) 



C ( b,N)C 2 (b, N) 
Cf(b,N) 

If N is large enough, and if we set b = b c (L, N = oo) we should get a value 
close to the number Q(b c , oo). We define an approximation to b c (L, N = oo), 
b c (L, N), by the equation: 

Q(b c (L, N), N) = ^/ 2( y = = 0.364739936 (5.22) 

The number above is taken from the same quantity defined for d = 2. 
Exercise: Calculate 0(6 = b c , oo) in d = 2. 
Now the objective is to establish that indeed the limit 

b c (L) = lim b c (L,N) (5.23) 



exits. We are now in a position to find the critical physical size l c again and 
make sure we get the same result as in the opposite order of limits. For this 
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we repeat the above procedure for several values of L. Then, we first invert 
the function b c (L): 

b c {L c {b)) = b (5.24) 

Actually, the set of values b c (L) is discrete, since L only takes discrete values. 
The function L c (b) defined above is a continuous interpolation of this set of 
discrete values. We now want to see that the following limit exists: 

l c = lim L c {b)a{b) (5.25) 

We find the same l c as before, with the same characteristics. To test for 
large N universality we need to go to a regime of asymptotically large values 
of N. We do this again first, before taking the continuum limit, assuming 
that this reversal of limits is allowed. 

We already know that Ojv(y, b) will exhibit critical behavior at b = b c {L) 
and y = as N — ► oo. There, it will obey large N universality if we can show 
that there exists a non-universal normalization factor, M(b, N, L), smooth 
in b at b = b c , such that: 



lim^oo M{b, N , L)0 N ( y = [^) '^ r) ,b = b c (L) 



i 

'l 



/3Na 2 (L) 



<(£,<*) (5.26) 

J\f(b, N, L) is a normalization factor we have to include just in order to get a 
finite answer. a±(L) and a,2(L) are iV-independent numbers, non-universal 
parameters. 

The universal function comes from carrying out the above calculation in 
d = 2 and is given by: 

C(£,a) = / due~ u ~ au +« u (5.27) 



Exercise: Verify the above and find out what a±, a,2 are in d = 2. 
Exercise: Calculate the function 0(a) in d = 2. 

The essential ingredient that needs to be verified now are the critical 
exponents 1/2,3/4. 

The parameter a2(L,N) is obtained by first setting 



b = b c (L, N) 



1 + 



/3Na 2 (L,N) 



(5.28) 



where b c (L, N) has been defined above. Next we write the derivative of 0, 
with respect to a at a = 0, which is at the critical size, and set the result 
equal to the corresponding universal number in the large N limit. 
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dQ(b,N) 



a=0 a 2 (L,N)V3N db 




b=b c (L,N) 

0.0464609668 



(5.29) 



Exercise: Verify the above formula. 

^ would be close to maximal at b = b c ; hence varies relatively little 
as b stays close to b c . Since b c is not known to infinite accuracy the reduced 
sensitivity on the exact value of b c is an advantage which motivates this 
choice for defining a,2(L,N). Unlike b c (L,N), the definition of a,2(L,N) 
involves going into the large N critical regime around b c (L, oo) and non- 
standard powers of N come in. If our hypothesis is correct, a,2(L, N) has to 
approach a finite non-vanishing limit ci2(L, oo) = a,2(L). This tests for the 
exponent 1/2. 

One also expects that the limit be approached as a power series in ^= . 
This is borne out by the data. 

To test for the other exponent, we first set 



and then form a ratio whose value at infinite ./V is again a universal number 
we can easily compute. 



This relation defines ai(L,N). 

Exercise: Verify the above formula. 

Similarly to 02 (L, N), the definition of a\(L, N) involves going into the 
large N critical regime around b c (L, 00). Consequently, we expect a\(L,N) 
to also have a finite limit a±(L, 00) = a\(L). This tests for the exponent 
3/4. The approach to the limit is again by a power series in 

Finally, to make sure that all this survives the continuum limit one needs 
to take L — > 00 at fixed physical size I = La(b c (L)), as explained earlier. 



A complete numerical test has been carried out in d = 3. The results 
are consistent with the large iV universality of the Durhuus-Olesen large iV 
phase transition in three Euclidean dimensions. However, it has turned out 
to be too difficult to determine the large TV" critical exponents from the data. 




(5.30) 




(5.31) 



5.3. Summary of numerical work 
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If we assume they are 1/2,3/4 we get consistency, but this is a somewhat 
weaker test then getting the exponents from the data at the expected values, 
with small errors. This type of numerical problem is common. 

In practice one does not need to use the definition of a(b) from the 
string tension as presented in the test. In d = 3 dimensions, as b — > oo, the 
asymptotic perturbative of a(b) is very simple: 

a(b) ~ c/b for b -> oo (5.32) 

There is a procedure to replace the constant c above by a function c(b) 
with c(oo) = c. The function c(b) is easy to extract from the Monte Carlo 
simulation, as it is directly related to the action. The large b asymptotic 
regime is entered for values of b that are much smaller than when c(b) is 
replaced by the constant c. This technique is called "tadpole improvement" 
or "mean field improvement". To test for the existence of a continuum 
limit, it is sufficient to use a(b) as defined above, and one does not need to 
also find S(6). One can find in the literature various determinations of S(6) 
which can be extrapolated to the range one uses and convert any continuum 
numbers that admit a dimensional interpretation into MeV units. There is 
no meaning to the actual number, as the string tension value we can assign 
a dimensional value to is the four dimensional one at N = 3. 

In d = 4, so far we have only a partial test: We know that smeared 
loops have a gap opening transition at a physical length at infinite N, but 
we have not tested for the critical exponents 1/2 and 3/4. Our guess is that 
consistent results similar to those in tree dimensions would be obtained. 



6. The bigger picture 

On the lattice, one can define another "string tension", which is loop 
size dependent, by using a Creutz ratio: 

ZcreuUL, b) = - log { {trW{LiL)){trW{L+hL + 1)) ) (6-1) 

where Wf(Li, L<i) is the Wilson loop matrix for an L\ x L<i- 

^Creutz{L,b~ l ) can be expanded for b^ 1 — > (the regular Feynman 
expansion) and for b — > (the so called "lattice strong coupling expansion" . 
The regimes of validity of these two expansions are disjoint; in between there 
is a crossover regime and we can bridge it only by numerical calculation. 

There are extra complications around 6 = 0. Rectangular loops of the 
type usually used have a "roughening" non analyticity in b at a point 6 roU gh- 
This non- analyticity is a lattice artifact. It can be avoided by choosing loops 
at generic angles with lattice planes. Then, the definition of EcVeutz needs 
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to extended. All this will increase the complexity of the strong coupling 
expansion. In the end, only a physical crossover separating the ranges of 
the weak and strong coupling expansion remains. We have no non-numerical 
calculational method to bridge it. To get the continuum string tension in 
units of the perturbative scale A.su(n) we need to take the continuum limit, 
which is, as we have seen, a correlated limit in which — > and the overall 
lattice scale of the loop goes to infinity. This correlated limit preserves the 
crossover. 

The idea we are pursuing is to improve the above scheme in two re- 
spects. First, since we wish to set up a calculation in continuum we forget 
about the lattice. Instead of thinking about ScVeutz we consider some other 
observable, for definiteness the extremal eigenvalue 9m of a Wilson loop of 
size I. 

For this to make sense, we need to be able to define 9 m in renormalized 
continuum field theory. We hope that this can be done by first constructing 
a renormalized polynomial in z corresponding to (det(z-Wf)) (for example 
by using smearing) and taking the roots of it to define 9m- While we have 
some idea how a calculation for small loops might proceed, for large loops 
we need something beyond ordinary field theory. Here we assume that an 
effective string model will describe (det(z — Wf)). This model will have a 
dimensional parameter, the string tension, and will be a good description 
for very large loops, with corrections parametrized by an increasing number 
of dimensionless parameters becoming more and more important as the loop 
shrinks. 

To relate the string tension to h-su(N)-, the dimensional parameter en- 
tering the perturbation theory for small loops, one needs to join the two 
regimes over the crossover. Here is the point that the simplification of large 
TV enters: At infinite N the crossover for zm collapses into a point and we 
have a phase transition. We postulate that we know that the transition 
is universal and that we know it is in the same universality class at the 
Durhuus Olesen transition. This postulate has good reasons to be correct, 
as I described. 

Therefore, for TV >> 1, the dependence of zm on intermediate scales, 
that is scales in the vicinity of the critical scale is known up to a few 
constants. This is the ingredient that was missing in the lattice scenario 
described above. It is now possible to imagine calculating to some order 
at short, intermediate and long scales and sewing together the three scale 
ranges. Requiring smooth matches could produce a number for the string 
tension in units of the perturbative scale A SU ( N y 

There are many variations possible, zm is only one possible example of 
a potentially useful variable, zm depends on the dilation of a fixed shaped 
Wilson loop, measured by I. As a function of I, zm will trace out a trajectory 
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from 9 = at I = (one could replace this by a < lA-su(N) <C 1) to = 
7r(l — 1/(27V)) at / = oo. For small I, the perturbative scale Agjy^) enters 
the calculations, and for large / the string tension enters. The two regimes 
are joined by the crossover. These lectures focused on this crossover, and 
this section has sketched the intuitive motivation for studying this crossover. 
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